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Abstract 

We consider the time-independent Wigner functions of phase-space 
quantum mechanics (a.k.a. deformation quantization) for a Morse po- 
tential. First, we find them by solving the ^-eigenvalue equations, using 
a method that can be applied to potentials that are polynomial in an 
exponential. A Mellin transform converts the ^-eigenvalue equations to 
difference equations, and factorized solutions are found directly for all 
values of the parameters. The symbols of both diagonal and off-diagonal 
density operator elements in the energy basis are found this way. The 
Wigner transforms of the density matrices built from the known wave 
functions are then shown to confirm the solutions. 



PACS: 03.65.-w, 03.65.Db, 03.65. Sq, 03.65.Nk 
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1 Introduction 



Phase-space quantum mechanics (see [2 [2], e.g.) is also known as deformation 
quantization (see [3]). In it, the density operator is replaced by its symbol, the 
Wigner function (or distribution) on phase space, and its equations of motion 
involve the Moyal *-product. Although a completely autonomous formulation 
of quantum mechanics results, the ^-eigenvalue equations obtained are difficult 
to solve. Few closed-form solutions have been found; see [TJ H] an d references 
therein. Here we derive a new, explicit solution of the *-eigenvalue equations of 
phase-space quantum mechanics by treating the Morse potential. 

The Morse potential is used in various physical applications. Our motivation, 
however, came from the deformation quantization of contact interactions. In 
particular, reference [5] uses the Morse potential to produce Robin boundary 
conditions at an infinite potential wall, one of the simplest contact interactions, 
as a limit of a smooth potential. 

In this paper we present the mathematical aspects leading to a "pure" defor- 
mation quantization of the Morse potential (see eqn. (fT2"]) below) with general 
coefficients. 

The exponentials in the Morse potential allow the pseudo-differential form 
of the *-eigenvalue equations to be replaced by a difference-differential equation. 
(Incidentally, the same is true for any potential that is polynomial in the ex- 
ponential function.) The difference-differential equations are then transformed 
into difference equations that can be solved for all values of the parameters. The 
properties of the Mellin and the inverse Mellin transform to relate derivatives 
with finite differences are extensively used. Details are given in the next section, 
where we also re-derive the known solution for the special case of a Liouville 
potential [3]. 

In the third section we provide an explicit check to show that the solution 
provided is indeed consistent with the Wigner transform of the density matrix 
of the Morse wave functions. Our final section is a short conclusion. 

2 Deformation quantization with a Morse 
potential 

Let us start with a brief review of Wigner functions and the Wigner- Weyl cor- 
respondence (see [H [2] , e.g.). For simplicity, we will restrict to one coordinate x 
and one conjugate momentum p, describing a flat two-dimensional phase space 
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R 2 . The generalization to M. 2N is straightforward. 

The Wigner function is related to the density operator of canonical quan- 
tization. More generally, every operator Q has a Weyl symbol Q(x,p) defined 

by 



-i(£x+r]p) 



Q(x,p) = WQ = — ^ / d^TrfQe" 

dy e~ iyp (x + hy/2\Q\x-hy/2) . (1) 



1 

2^ 

The map W : Q i— > Q{x,p) 1 from operators to phase space functions (and 
distributions) is called the Wigner transform. It is a homomorphic map from 
the algebra of operators to the *-algebra of symbols: 

W{QR) = W{Q) * W{R) . (2) 

Here the symbols are multiplied using the Moyal *-product, 

* = exp | y { d x d p - d p d x ^ | (3) 

for consistency with the Wigner transform (TTJ). An inverse W _1 also exists — it 
is commonly referred to as the Weyl map: 

Q =: W~ l Q{x,p) = — — r / d£ dr] dx dpQ(x,py^ £ - x)+l ^P- p K (4) 

Up to normalization, the Wigner function p(x,p) is defined as the Wigner trans- 
form of the density operator p: 

Pi*,P) ■■= ^ W (P) ■ ( 5 ) 

The Wigner function is a quasi-probability distribution. Expectation values 
are calculated as 

(Q) = J dxdp p(x,p)Q(x 7 p) . (6) 
However, even with a pure state density operator p = \ip)(ip\, 

p^p) = 1 d y e ' typ ^ x + h y/^ W x - h yl 2 ) ( ? ) 

{2-Kyn j_ OQ 

is negative somewhere (except for the special case of Gaussian wave functions) . 
The Wigner function satisfies the equation of motion 

ihd t p(x,p,t) = [H,p(x,p,t)] r , (8) 
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where [H, p]^=H*p — p*H. It can be expressed as a linear combination of 
stationary Wigner functions with time-dependent coefficients: 

p(x,p,t) = £ C ELER e-^- E «Wp ELER (x,p) . (9) 

Here P E E = W(\El)(E r \)/2tt?i denotes the Wigner transform of a matrix 
element of the density operator in the energy basis. As *-eigenfunctions, they 
can be found by solving the system of equations: 

H *Pe l e r ( x 'P) = El Pe l eS x 'P) ' 

P ElEr (X'P) *H = E r P ElEr (x, P ) . (10) 
These are known as the ^-eigenvalue (or sometimes "stargenvalue" ) equations. 
Alternatively, the Wigner transform 

/oo 
dy (x + %y/2\E L ) (E R \x - hy/2) (11) 
-oo 

allows them to be determined from the wave functions, if known. In the case of 
smooth potentials, the resulting Wigner functions are known to agreed 

The goal here is to perform a "pure" deformation quantization by solving the 
^-eigenvalue equations directly, without reference to operators or wave functions. 
This will be done for the Morse potential 

?>V 

V(x) = ( e~ 2ax - e~ ax ) . (12) 

2m 

With its short range repulsion and longer range attraction, this smooth po- 
tential has been useful in many physical applications. In particular the Morse 
potential can be used to recover infinite wall with Robin boundary conditions 
in deformation quantization [5]. We will use the Mellin transform to convert 
the *-eigenvalue equations to difference equations. 

We use a new method that produces difference equations for potentials that 
are polynomials of an exponential in x0 This is significant since solutions to 
the ^-eigenvalue equations are generally difficult to find (see [4], e.g.). To see 
why, let us start with a more general Hamiltonian and then specialize to the 
Morse potential. Writing 

El —'■ ~z i E R =: ; 13 

2m 2m 



For discontinuous potentials, however, that is not necessarily the case 6 7 8 5. In [5] 
the results reported here are applied to study the example of the infinite wall, or equivalently, 
a particle confined to the half-line. 

2 Strictly speaking, a difference equation is obtained for any potential that is a linear 
combination of exponentials, exp(— ctix), i = 1, . . . ,n, say. That is not likely to be helpful, 
however, unless all the ratios ai/aj are rational. 
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the ^-eigenvalue equations (|10[) arc of infinite order in momentum-derivatives 
for a generic Hamiltonian H — p 2 /2m + V(x): 

/ n=0 v 7 

and similarly for the "right" equation. However, the exponential form of the 
Morse potential (fT2j) allows them to be written as differential-difference equations- 
the exponentials generate translations in the momentum. Their explicit form 
becomes 

h 2 - ihp 

^ 2 xP+ -0 xP = (15) 
&V _ 2ax . . /3ftV _ Q;E / iha\ , f P 2 -h 2 k 2 \ 

and 

8^»-P-2^^= ( 16 ) 

ftV _ 2q:c . . p% 2 K 2 / iha\ fp 2 -h 2 k 2 R \ 

— e ^p( x , P + i na)- P —e^ P ( X , P+ —j + (^^y. 

The integral transform technique leads to further simplifications. Suppose 
the Wigner function can be written as 

p(x,p) = R(u,p) , u := We^a^/n 4 . (17) 

The Mellin transform of the Wigner function is 

/>oo 

W{s,p) := M{R}(s 7 p) = / u 3 - 1 R{u lP )du . (18) 

Jo 

To transform (fT5j) and (fT6|) into difference equations for W(s,p) we consider the 
inverse Mellin transform 

i />c+ioo 

R(u,p) = — / u- s W(s,p)d S , (19) 

where the constant c can be any constant for which the transform converges 
according to the Mellin inversion theorem. 



2.1 Solution for the Liouville potential 

We will first apply the method to the Liouville potential 

V L {x) = t± e ^ . (20) 
2m 
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It is the simplest case as it can be viewed as a Morse potential for b = 0. Also, we 
can check our results since the Wigner functions for this potential have already 
been found in 0]. 

The ^-eigenvalue equations (|10p imply the difference equations 

(p/h + 2ias) 2 W (s,p) + (2a) 2 W (s- l/2,p-iah) = k 2 L W (s,p) , 
(p/h-2ias) 2 W a (s 1 p) + (2a) 2 W a (s- 1/2, p + iah) = k R W (s,p) . (21) 

Let us now assume that the solution is factorized into two parts - "left" and 
"right" factors - each depending on fcz, or k R only: 

W(s,p) = Nw L (s-^,k L )w R (s+^,k R ) , (22) 

with N a normalization constant. For the left factor we find 

(p/h + 2ias) 2 w L (s - ip/2ah) + (2a) 2 w L (s - ip/2ah - 1) = (23) 

= k\ wl(s — ip/2aK) 

and for the right factor 

(p/h - 2ias) 2 w R (s + ip/2ah) + (2a) 2 w R (s + ip/2ah - 1) = (24) 

= k 2 R w R (s + ip/2ah) . 

Using the substitution t = s — ip/2ah we arrive at 



w(t - 1, k L ) 



w(t,k L ) . (25) 



(2a) 2 _ 

Equation (|24p also leads to the above equation if we use t = s + ip/2ah and 
k R instead. Therefore we need only work with (j2"5j) and the solutions will just 
differ in their arguments and labels of k. 

The solution of (123 is 



w(t, k L ) = T(-t + ik L /2a) T(-t - ik L /2a) , (26) 

by the defining property T(z + 1) = ^r(z) of the gamma function. Tracing back 
to equations (I22H23[) . (fT!?)) and p7|) . we can write the Wigner function in terms 
of the inverse Mellin transform: 

r-C+lOQ 

P khkR {x,P) « / ds u- s x (27) 



c — ioc 



i(p/h±k L )\^( i( P /h±'k R ) 

-s 



, 2a I \ 2a 



G 



This last is an integral representation of the Mcijer G-function. Using equation 
(43) on pg. 353 in [S] the Wigner function becomes 



k L k_ 



(28) 



G 



i{p/h + k L ) i(p/h~k L ) i(p/h-k L ) i(p/h + k L ) 



2a 



2a 



2a 



2a 



where we used the identity 



G 04 (u\l - ai, 1 - a 2 , 1 - a 3 , 1 - a 4 ) = G 04 (l/u|oi, a 2 , a 3 , a 4 ) 



(29) 



As it should, the formula for the Wigner function (|2T)1) coincides with the one 
obtained using different methods in [4 . It describes the phase-space quasi- 
distribution for a particle in a Liouville potential. The advantage of the method 
proposed here is that it can be generalized to the Morse potential. 



2.2 Solution for the Morse potential 

Now let us go back to the original problem of finding the Wigner function for 
the potential (fT2"j) . The left ^-eigenvalue equation has the form 

(p/h + 2ias) 2 W b (s,p) + (2a) 2 W b (s - 1/2, p - iah) (30) 
-^(2a) 2 W b (s- l/4,p-iah/2) = k 2 L W b (s,p) , 

where we set b = /3k/ 'a . We also have the complex conjugate (right) equation, 
with fci replaced by fc^. To solve the new left difference equation (|30|) , we 
substitute the ansatz (|2"2"]l to obtain 

(p/h + 2ias) 2 w L (s - ip/2ah, k L ) + (2a) 2 w L (s - ip/2ah - 1, k L ) (31) 
--(2a) 2 w L (s - ip/2ah - 1/2, k L ) = k 2 L w L (s - ip/2afi,k L ) . 

Using the same substitution t = s — ip/2ah as in (|24l) and switching to w b = wl 
to account for the parameter dependence, we arrive at the difference equation 
relevant to the Morse potential: 

w b (t-l,k L )- ~w b (t -l/2,k L ) = 

This equation has a trivial solution for 6 = 0, the Liouville case. The right 
factor satisfies an identical equation with replaced by &r and t = s + ip/2ah. 



t 2 + ■ kh 



(2af 



w b (t, ki 



(32) 
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For 6 = 1 the solution can be written in terms of gamma functions, as in the 
6 = case 



ir, .(/./,-) x r | -t + Y) r ( -* ' 



2a 



(33) 



1 ik , ^ , 

-t h — i — r -t- 

2 2a) V 2a 



The inverse Mellin transform then gives us the Wigner function: 
Pi 



,-i40 
^04 



Lr Q4 



,-f40 
^04 



(a;,p) oc 


(= 


ip 




2aE 




ip 


(= 


2aE 




1 


(= 


2 + 




1 


(= 


2 + 



(34) 



2a' 2 
ife 1 
2a '2 



ip 

2ah 
ip 
2ah 



ip 

2aE 

ip 



ikL ip 
2a ' 2oft 



ifct 
2o~' 

2o~' 

ifcz, 
2a ' 



ip 

2oft 



ip 

2ah 



1 



ikn 1 
2a' 2 " 
ip ifcfl 
2oft + 2a '' 
ikn 1 _ 
2a ' 2 " 
ip i/c_R 



u 2 2a?i 2a ' 2a?i 2a ' 2 2a?i 2a ' 2a?i 2a 



ip 


ikR 


2ah ~ 


2a 


ip 


ikR 


2ah ~ 


2a 


ip 


ikR 


2ah ~ 


2a 


ip 


ikR 


2aE ~ 


2a 



Another solution that is easy to find is for 6 = 2:|f| 

,(*, k) « (t + 1/4) nr (-t ± I) - n r (-* + 5 ± | 



(35) 



It can be written in a different form which shows a pattern shared with the case 
6 = 1: 

w 2 {t, k) cx (2ifc/a + 1) T (-t + ik/2a) T (—t + 1 — ifc/2a) + (36) 
4ifc 

r (-t + 1/2 + ifc/2a) r (-t + 1/2 - ifc/2a) + 

a 

(2ifc/a - 1) T (-t + 1 + ifc/2a) T {-t - i/c/2a) . 

The Wigner function can be found with a trivial but lengthy calculation that 
is essentially identical to the 6 = and 6 = 1 cases. The exact combination of 
Meijer G-functions is not of interest to us; we will derive a general expression 
that includes this one later. 



A useful observation is that the left factors (for different 6) can be written 



as: 



wi = Ciw (t- l/4,k + ia/2) + C 2 w a {t- l/4,fc-«a/2) 



(37) 



3 The 6=2 Wigner function can also be found from the Liouville case using supersymmet- 
ric quantum mechanics. In deformation quantization, the ladder operators of supersymmetric 
quantum mechanics are replaced by functions and a star product is used, however the transi- 
tion is fairly straightforward [I] . 
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and 

w 2 = C x w {t- 1/2, k + ia) + C 2 w (t- l/2,k) + C 3 w (t- 1/2, k-ia) . (38) 

This suggests that by choosing the constants correctly we can write the solution 
for any integer b as 

w b (t,k)oc ^2 C*w (t-b/4,k-ina). (39) 
n=-| ,-§+i,...,§ 

We can substitute this ansatz into the equation ()32j) using undetermined coef- 
ficients. In principle, comparison of the coefficients of independent terms can 
determine C\ for any b. This seems to fail, however, in the case of non-integer b. 
Furthermore, even for the simplest cases this program is very difficult to carry 
outH Clearly we need an algorithm that reproduces the constants directly and 
allows a generalization to include all Morse potentials of the form (fT"2"|) . 



2.3 Systematic solution of the difference equations 

We now show how to find the relevant solutions of the *-eigenvalue equations for 
the Morse potential, for all b € R+. We exploit once again the property of the 
Mellin transform to relate differential equations and their solutions to difference 
equations and their solutions. 

To convert our difference equation (|32p into a differential equation we use 
the following two properties of the Mellin transform: 

M{r 2 f" (r) + Tf(r)}(s) = s 2 M{f(r)}( S ) , (40) 
M{T a f(r)}( S )=M{f(r)}( S + a). (41) 

If we apply the inverse Mellin transform directly to (|32p . we end up with an equa- 
tion that we cannot solve. This is because the Mellin transform converts argu- 
ment translations into powers of the argument via (|41|) . To eliminate fractional 
powers, we use the substitution s = 2t. The new equation for Wb(s) = w(t(s)) 



w b (s - 2) - ^w b (s - 1) 



2 



21 {2af 



Ms) (42) 



results in a simpler, integrable equation: 

r 2 f"{r) + rf{r) + [{k/af - 4/r 2 + 2b/r]f(r) , (43) 
where w(s) = M{f(T)}(s). 



4 For integer b, see eqn. I|52| l below, however. 
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The solution /(r) of this equation can be found if we make the substitution 
/(t) = T~ 1 / 2 g(r) and then u(z) = g(t(z)), where z = 1/t. The new function 
u(z) satisfies the so-called Whittaker equation, treated in [TU], Chapter XVI, 
and also in the Appendix. Its two linearly independent solutions are defined 
in [TT], pg. 755. They are called Whittaker functions and can be expressed 
in terms of the Tricomi confluent hypergeometric function U(fx, v 1 z) and the 
Kummcr confluent hypergeometric function M(/i, 

M lm (z) = z ™+i/2 e -^/2 m (1/2 + m - l t 1 + 2m; z) , (44) 
Wi m {z) = z ™+ 1 /2 e - 2 /2 17(1/2 + m - 1, 1 + 2m; z) . (45) 

For our purposes, we only need the definitions of those functions 

M (/W) = £^, (46) 
- 1 



n=0 v ' 



U(n, v; z) = r( ^, ^ z>- v M(\ + n-v,2-v;z) 



r " " ] -M{n,v;z) . (47) 



rOu-^ + i) 

Here we use the Pochhammer symbol (fJ.) n '■= A*(m + + n — 1), (/i)o := 1. 

The solution is: 

/(t) =C , ir 1/2 Mi^(4/r) + C , 2 r 1/2 ^| m(4/t) . (48) 

This is the general solution and it therefore depends on two arbitrary constants, 
Ci, C*2- We must set C% = to describe the physical states, however. To see 
that we have to transform back to the solutions of the difference equations and 
compare with the known solutions, re-derived in the previous sections. 

Let us first confirm that (|48|) . with C\ — 0, indeed recovers the known 
solutions (f2l>l 1551 156]) for b = 0, 1, 2, respectively, and that it also justifies the 
ansatz (|55|) for all non-negative integer b. For 6 e No we can write the Whittaker 
VF-function in terms of the modified Bcsscl functions 

y 2 I Y — n 



w±M = —F=r hr + H x ( 49 ) 



" (_l)»+fc( 2 fc - n + 2 M )(-n)„_ fc 
X ^- o r(n-A;)(fc-n + 2 A i)„ +1 ! - 



5 The Whittaker function M[ m (z) should not be confused with the Kummer function 
M(/i, v, z) in the above equation. Subscripts are used to denote the parameters of the Whit- 
taker functions in the literature, and the explicit bracket notation is used for confluent hyper- 
geometric functions. For further information involving the hypergeometric functions see 1111 . 
p.753 and [12], p.503-506. 
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With the help of the integral representation of the Bessel functions (effectively 
finding the inverse Mellin transform), 

K V{Z ) = -I -„)(-) d s, (50) 



we can find the solution of (|42|) for integer b: 
b 



w b 



(t, fc) « E C h n T (-t + n/2 + ik/2a) T (-t - n/2 + 6/2 - ik/2a) . (51) 



n=0 



This is nothing more than (|39p with a shifted summation index. The coefficients 
are now explicit, however: 

r b = (-ir(2n-b + 2ik/a)(-b)^ n 

(b-n)l(n-b + 2ik/a) b+1 ' 1 J 

The Wigner function for b G No can be found from (|5"Tj) as in the preceding 
two subsections, with result 

b 

PkLkR (x,p)= J2 C b m C b n x (53) 



m,n— 



xGq4 f - 



n i(p/h + kL) b n i(j>/Ti — ki/) 
2 + 2^ '2 ~ 2 + 2^ ' 

m i(jp/fi — kft) b m i(p/h + kfj) 
~2 2^ '2 ~ 2" 2^ 

Bound states can also be treated this way; one only needs to consider imag- 
inary k, for energies E < 0. The resulting form of the Wigner function differs, 
however, from the expression found by the Wigner transform 

(54) 

To write it in this form we use the properties of the Whittakcr functions. Recall 
that the energies are given by (|76[) and we can write 

/(r) cx T 1/2 Wt ^(4/t) = t 1/2 W^^_ u _i{A/t). (55) 

The relationship 

V-iW = (-!)"«/! ^-"e-*/ 2 ^- 2 "- 1 (56) 
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(for integer v) and (fT5j) allow us to write the solution as 

Transforming the above we find the corresponding factor: 

- (— 9) l 9 2t (b — v — 1\ 
w(f)oc^LJ _ )r(-2t + Z + 6/2-i/-l/2). (58) 



With the help of the inverse Mellin transform we find ([5"4")l as in the unbound 
case, using (f50|) and ([19]). 

Let us proceed to the case of non-integer b. For b € K, we find an explicit 
and closed expression for the solution of the difference equation. We start by 
rewriting the solution of the differential equation (|43j) in terms of hypergeometric 
functions: 

f{r) = e-^r-^u(\- | + **,!+ «*; 1) . (59) 

This is necessary in order to perform the inverse Mellin transform of /(r) in 
closed terms, which presents a technical problem if we use the Whittaker func- 
tion. 

We use the relationship (1471) between the Kummer and Tricomi hypergeo- 
metric function and the integral expression 

/>C+lOO 

e~' Ta: M(j3,r,Xx)= ds x~ s a~ s r(s) 2 F 1 (f3, s ;j ; Act" 1 ) (60) 

J c— zoo 

to make it possible to find the inverse Mellin transform of (|5"9")l . It is given in 
term of the Gauss hypergeometric function 2-Fi- Switching back to our original 
variable t, and indicating the fc-dependence explicitly, we can write 

4 t+lfe / 2Q r(-2iA:/a) n , . , 
^ fc) - r(l/2-V 2 -^/a) r( ^ + ^ /Q)X 

2-Fi (1/2 - 6/2 + ifc/a, -2i + ik/a; 1 + 2 ife/a; 2) + 

4t-ik/2* T (2ik/ a ) 



-T(-2t - ifc/a) x 



T(l/2-b/2 + ik/a) 

2 F 1 (l/2-b/2-ik/a,-2t-ik/a;l-2ik/a;2) . (61) 

This is the solution of the difference equation (|32|) for any real b, including the 
ones we already found. After substitution in (|ITI[l"9ll2"2"j) . it yields our main 
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result: 



l_ ioa dsu ~ Swb { s -^k M ) Wb { s+ ^ kR ) ■ (62) 



3 Wigner functions from wave functions for a 
Morse potential 

For completeness, let us confirm that our solutions of the *-eigenvalue equations 
(fT0|) are the same as the Wigner functions derived from the wave functions (see 
the Appendix) using (fTT]) . 

Consider the unbound states first 

m = Ce-y/z Ay^« M (I - b - + ^,l + ^;y) + 

\2 2 a a J 

+ Ce-y/ 2 A*y- ik /«M(]--^-^,l-—iy) ■ (63) 
\2 2 a a ) 

Using the substitutions w = e~ ahy / 2 and v = 2ne~ ax /a we can rewrite the 
integral transform The result involves integration over w of products of 

the type 

exp[-v(w + w- 1 )/2]w m v n M (a^b^vw) M (a 2 ,b 2 ;v/w) , (64) 

which (to the best of our knowledge) are not integrable in closed form. However 
we can expand the Kummer M -function as in (|46[) . Then all the integrations 
can be performed explicitly using the integral representation of the Bessel K- 
function 

K u {z) = - dw w -(*+i) e -i*(v>+iM . (65) 
2 Jo 

This procedure leads to an infinite sum: 

P (v,p)oc V ^ [ iV* fc /» { x\ m{ ^ n K n _ m _ 2ip/na (v) + (66) 

^ m, "< (?)m(c)r 



m\n\ 

m,n— 



I 7|2 \X)m\X)n v / \ . 

\A\ y-r — p— J^-m-n+2i(k-p/h)/a( v ) + 

I Af2 [X)m{x)n js r \ . 

\A\ -pr —y ft m -n-2i(k+p/h)/a(V) + 

\S )m\^ )n 

I A*\2„-2ik/a (X)m(x)n v / \ 

\S)m(S)n 
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with x '■= 1/2 — b/2 + ik/a and <; := 1 + 2ik/a. This expression can now be 
used to compare with (|61l [62]) . 

To do that we need to calculate explicitly the contour integral (|62|) . We are 
facing a similar problem — the integrand is too complicated and we need to 
rewrite the Gauss hypergeometric function as 



in order to integrate in closed form. With the use of the contour representation 
of the Bessel if-function and some simple algebra we recover the infinite 
sum This confirms that the Wigner function (|6"2")l indeed coincides with 

the Wigner transform of the density matrix for the Morse potential from the 
Schrodinger treatment. 

The Wigner functions of the bound states can be obtained from the known 
wave functions (|74|) using the integral transform (fTTjl FI Equation (|75|) allows 
us to evaluate the integral in closed form using the modified Bessel functions 
K v (x). The substitutions v — 2ne~ ax /a and w — exp(— othy/2) in (fTTj) and the 
integral representation (|55t are used. 

4 Conclusion 

Our main result is the solution (|62l [B"Tj) of the ^-eigenvalue equations (flOl) for 
the Morse potential (fT2j) with arbitrary real b. It subsumes the simpler formula 
([5511551152]) valid for all b e N . 

The solutions obtained have already been applied to a study of Robin bound- 
ary conditions in phase-space quantum mechanics [5j. 

It should also be possible to use our method to solve the *-eigenvalue equa- 
tions for other potentials that are polynomial in an exponential (say exp(— ax)). 



6 For 6 = 2, the bound-state Wigner functions for the Morse potential il 1 21) have already 
been treated this way in | 14| . However, they are not obtained there by solving their dynamical 
equations, as we have done. In addition, while we find the unbound-state Wigner functions 
in the same way, the unbound states are not considered in [14] , 




(67) 
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Appendix: The Morse potential in Schrodinger 
quantum mechanics 



Following Matsumoto |13] , we can solve the stationary Schrodinger equation 
for the unbound wave functions of the Morse potential ([T2"]l . The substitution 
?p(x) = (f>(z), z = exp(— ax), changes the Schrodinger equation into 



zV + zcf>' + 



1 



2mE 



K 2 Z 2 + K Z (i Z 



cj> = 



(68) 



This can be further transformed into canonical form (without a first derivative 
term) using the substitution 0(z) = z -1 / 2 F(z). Changing the variables to 
y := 2nz/a leads to the Whittaker equation, treated in [ID] . Chapter XVI: 



/* 



b 1 

2y 



IT 



f = , 



(69) 



where f(y) := F(ay/2n) and, as before, k — y/2mE/% and b — /3n/a. 
Now the wave function can be written as 



c/2 



CiM k in{y(x)) + C 2 W>, ih.{y{x)) 



(70) 



Imposing reality yields C% — 0. The second term has physical asymptotic be- 
haviour: for large positive x it is sinusoidal with a phase depending on the 
potential parameters; for negative x far from the origin, there is the expected 
rapid exponential decay of a classically forbidden region. The wave function is 
therefore 

(71) 



With the help of equation (|47l) we can rewrite this result in a form similar 
to that given by Matsumoto in [13] for a Morse potential with 6 = 2. The wave 
function is manifestly real in this form: 



il,(y) = Ce-y/ 2 Ay^M ( I - b - + ±,1 + ™ 



\y 



Ce-y' 2 A* y- ik ' a m(-----,1-— ; , 
V 2 2 a a 



2 2 

with C a real normalization constant, and 

p/ 2ik \ 

A _ ^ a I 



ik \ 



(72) 



(73) 
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Let us now consider the bound states. Their wave functions are given in j!4| 

as 

i/>(x) oc exp(- Ke - Q " / a ) e -a(»-b/2+i/2)x L b - 2 "- x (2Ke- ax /a) , (74) 

where 

*<•>-£<-'>-(;:*)£ to 

m=0 v 7 

are the associated Laguerre polynomials L^(x). The energies are 

E v = -^-{ V -b/2 + l/2f, (76) 
for integer ^ £ [0,\_b/2\], where [a\ is the smallest integer less than a. 
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